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Abstract 

We use the Constitution supernova, the baryon acoustic oscillation, the cosmic microwave back- 
ground, and the Hubble parameter data to analyze the evolution property of dark energy. We 
obtain different results when we fit different baryon acoustic oscillation data combined with the 
Constitution supernova data to the Chevallier-Polarski-Linder model. We find that the difference 
stems from the different values of f2 m o- We also fit the observational data to the model independent 
piecewise constant parametrization. Four redshift bins with boundaries at z = 0.22, 0.53, 0.85 and 
1.8 were chosen for the piecewise constant parametrization of the equation of state parameter w(z) 
of dark energy. We find no significant evidence for evolving w(z). With the addition of the Hubble 
parameter, the constraint on the equation of state parameter at high redshift is improved by 70%. 
The marginalization of the nuisance parameter connected to the supernova distance modulus is 
discussed. 
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I. INTRODUCTION 



Since the discovery of the late time cosmic acceleration by the Type la supernova (Snla) 
observations jl, 2], a lot of efforts have been made to understand the driving force behind the 
cosmic acceleration. The standard models in cosmology and particle physics give no answer 
to this problem. To address the problem, one needs to modify either the left hand side or 
the right hand side of Einstein equation. Modifying the left hand side means that general 
relativity is modified, models such as the Dvali-Gabadadze-Porrati model [3j, and f(R) 
gravity [4|-|l3| , have been proposed along this line of reasoning. On the other hand, in the 
framework of Einstein gravity, an exotic form of matter with negative pressure, dubbed as 
dark energy, has to be introduced into the right hand side of Einstein equation to explain the 
phenomenon of cosmic acceleration. However, the nature and origin of dark energy remain a 
mystery. Many parametric and nonparametric model- independent methods were proposed 



to study the property of dark energy, see for example 



14M34J] and references therein. 



Recently, it was claimed that the flat ACDM model is inconsistent with the current data 



at more than la level 



35M38| . Furthermore, it was su gges ted that the cosmic acceleration 



is slowing down from z ~ 0.3 in 35] . The analysis in 35] is based on the commonly used 
Chevallier-Polarski-Linder (CPL) model 33|, |34| with the Constitution Snla 39] and the 
baryon acoustic oscillation (BAO) distance ratio data 40], the result is somewhat model 
dependent. In [36], the authors applied the Union Snla data 41], together with the BAO 
A parameter 42] and gamma-ray bursts data to the piecewise constant parametrization 
of the equation of state parameter of dark energy. In this analysis, the equation of state 
parameter w of dark energy is a constant in a redshift bin, and three redshift bins with 
boundaries at z — 0.2, 0.5 and 1.0 were chosen. They also analyzed the Constitution Snla 
and the BAO distance ratio data. In 37], the authors considered two redshift bins by 



using the Constitution Snla data, and t hey found that dark energy suddenly emerged at 
redshift z ~ 0.3. The results obtained in [38] were based on the analysis of the Constitution 
Snla, the full Wilkinson microwave anisotropy probe 5 year (WMAP5) (43! and the Sloan 
digital sky survey (SDSS) data. For the redshift range < z < 1, four evenly spaced 
redshift bins were chosen. However, choosing the same four evenly spaced redshift bins, the 



authors in 



44J found no evidence for dark energy dynamics by using the Constitution Snla, 



WMAP5, BAO (45], the integrated Sachs- Wolfe effect, galaxy clustering and weak lensing 
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data. Different data sets and analysis may give different results. In this paper, we apply the 
piecewise constant parametrization of the equation of state parameter w(z) to do a more 
careful model independent analysis. We choose four redshift bins by requiring NAz ~ 30 



in each bin, and we use the Constitution Snla |39|, the BAO |40|, the derived WMAP5 43 1 



and the Hubble parameter H(z) data 46|-|48|. 



This paper is or gan ized as follows. In section II, we first review the analysis by using 
the CPL model in [35 1 , and find that their result heavily depends on the choice of BAO 
data. By using the BAO distance ratio data, the best fit value of fi m o — 0.45l ' 1:1 , which 
is not consistent with other observational result. On the other hand, the best fit value 
of fl m Q = 0.291qq4 if the BAO A parameter is used. Then we apply the Snla, BAO, 
WMAP5 and H(z) data to study the property of dark energy by using the piecewise constant 
parametrization of the equation of state parameter w(z), in section III. We conclude the 
paper in section IV. 



II. OBSERVATIONAL CONSTRAINTS ON CPL PARAMETRIZATION 



To study the dynamical property of dark energy by observational data, one usually pa- 
rameterizes the equation of state parameter w(z). Following [35|, we first study the CPL 
parametrization 



w(z) = Wq + 



WgZ 

l + z 



(1) 



The dimensionless Hubble parameter for a flat universe is 
H 2 (z) 



E\z) 



H 2 



a m0 (l + zf + (1 - n m0 )(l + zf(i+-o+^) exp(-3?/W(l + *)). (2) 



In this model, we have three parameters Q m o, wq and w a , let us denote them as p 



(£l m0 ,Wo } w a ). We first use the Constitution compilation of 397 Snla data [39[ to constrain 
the model parameters p. The Constitution sample adds 185 CfA3 Snla data to the Union 
sample 4l|. The addition of CfA3 sample increases the number of nearby Snla by a factor 
of roughly 2.6 — 2.9 and reduces the statistical uncertainties. The Union compilation has 



57 nearby Snla and 250 
the ESSENCE Survey [50J 



iigh-,2 Snla. It includes the Supernova Legacy Survey 49| and 



5l|, the older observed Snla data, and the extended data set of 
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distant Snla observed with the Hubble space telescope. To fit the Snla data, we define 

397 r / \ / rr n \l2 



XsnKPi h 0> = 2^ 2 

i=l Gi 



(3) 



where the extinction-corrected distance modulus \i(z) is the difference between the apparent 
magnitude m(z) and the absolute magnitude M of a supernova at redshift z, 

fi(z,p,H^) = m(z) — M = 25 - 51og 10 # n + 5 \og 10 [D L (z)/Mpc}, 

the absolute magnitude M applies equally to all magnitude measurement, and its effect is 
manifested by the nuisance parameter H$ ; a\ is the total uncertainty which includes the 
intrinsic uncertainty of 0.138 mag for each CfA3 Snla, the peculiar velocity uncertainty of 
400km/s, and the redshift uncertainty [391]; and the Hubble constant free luminosity distance 
D L (z) = H d L (z) is 

l + z . 



D L (z lP ) = H d L (z) 



=smn 



where 



sinn 



(VP 



sin(y/\Q k \x)/y/\^ 



■r. 



smh(y/\^\x) I \/\fk 



dx 
E(x,p) 

if Q k < 0, 

if n k = o, 

if tt k > 0. 



(4) 



(5) 



The nuisance parameter Hq is marginalized over with a flat prior when we apply the Snla 



data. For the details of the marginalization method, see 26 



23| ■ Then we add the BAO 



parameter A = 0.469(0.96/0.98) 0,35 ±0.017 42] into the Snla data to determine the param- 
eters p. So we have X 2 (p) = [A- 0.469(0.96/0.98)-°- 35 ] 2 /0.017 2 + xL- The BA0 parameter 
A in a spatially flat universe is defined as 



A(p) = y/ti 



H D v (z bao ,p,H ) VU 



niO 



mO 



Zbao %bao 

where z bao = 0.35, and the effective distance 



Zba 



E(z ba 



^bao 



dz 
Eh) 



1/3 



D v (z,p,H 



1 



D 2 L (z) 



1/3 



(6) 



(7) 



[l + z) 2 E(z 

Finally we add the shift parameter R with which the Z-space positions of the acoustic peaks 
in the angular power spectrum shift, to the combined Snla and BAO A data. The shift 
parameter 

(8) 



pi 

R(p) = A/^mO / 
JO 



E(z) 



1.710 ±0.019, 



where the last scattering surface redshift z\ s = 1090.0. So now we minimize 

(R - 1.71) 2 [A - 0.469(0.96/0.98)-°- 3512 



x 2 (p) 



0.019 2 



0.017 2 



+ v 2 . 

1 Asn 



Fitting the Constitution Snla data, the BAO distance ratio Dy(z = 0.35)/ Dy(z = 0.20) = 
1.736 ± 0.065 (hereafter BAO I) [45] and the shift parameter R to the CPL model, the 
authors in 35] conclude that the functional form of the CPL ansatz is unable to fit the 
data simultaneously at low and high redshifts. In this section, we replace the BAO distance 
ratio data (BAO I) by the BAO parameter A (hereafter BAO II) to fit the CPL model. By 
fitting the Snla, Snla+BAO I (II), and the Snla+BAO I (If)+R data to the CPL model, 
we obtain the joint constraints on the parameters Q m o, wq and w a in the CPL model. 
The values of Q m Q and x 2 constrained from different data sets are shown in Table I. For 
the purpose of comparing the results, we only show the constraint on Q m o in Table I. For 



the Snla+BAO I data, the joint la constraints are fl m o = 0.45_ ail , wq = — 0.13_ ; 95 and 
w a = — 12.2l}5 3. Fitting the Snla+BAO II data to the CPL model, the joint la constraints 
are fl m0 = 0.29^04, Wq = — 0. 901^37 an d w a = — 0.6*33 ■ Comparing these results, we find 
that Wq and w a are consistent with each other at la level, but Q m Q is barely consistent with 
each other at la level. For the Snla data or the Snla+BAO I data, the best fit value of 
fl m o is much larger than the value Q m o ~ 0.3 obtained from other observational constraint, 
and that makes the result inconsistent with the ACDM model at more than la level. The 
BAO II data or the WMAP5 data lowers the value of Q m n, and therefore makes the resu 



consistent with the ACDM model at la level. To compare our results with those in 35], 



Data 



Snla 

Snla + BAO I 
Snla + BAO II 
Snla + BAO I + R 
Snla + BAO 11 + R 



CPL model 



n 



m0 



0A5 +0 il X 2 = 462.07 



ft m0 = 0.45i°;?I, x 2 = 462.44 
n m0 = 0.2918:81, x 2 = 466.18 

n m0 = o.26i8:8l» x 2 = 467 - 74 

n m0 = 0.27 ± 0.03, x 2 = 466.81 



ACDM model 





= 0.29 + 0.02, 


x 2 


= 466.32 


^mO 


= 0.29 + 0.02, 


x 2 


= 467.61 


^mO 


= 0.28 + 0.02, 


x 2 


= 466.42 


^mO 


= 0.27 + 0.02, 


x 2 


= 469.49 


^mO 


= 0.27 + 0.01, 


x 2 


= 468.42 



TABLE I: The la error estimate of O m o for the CPL model and ACDM model, 
we also show the evolutions of q(z) and Om(z) constrained from Snla, Snla+BAO II and 



Snla+BAO ll+R in Fig. [TJ The deceleration parameter 

. _ n mQ (l + zf + [1 + 3w(2)](l - Ono)(l + *) 3 < 1+w ° +w <0 exp(-3w a z/(l + z)) 
q{z) 2E^) ' (9) 

and Om(z) is defined as 

E 2 (z) - 1 

- (TT^r (10) 

The evolution of q(z) gives us the information about how fast the Universe expands. The 
sign of q(z) shows whether the expansion is accelerating or decelerating, q(z) > means 
deceleration. For a ACDM model, Om(z) = Q m o is a constant. Also at low redshift, larger 
Om(z) means larger w of dark energy. 

C omp aring the evolutions of q(z) and Om(z) shown in our Fig. [I] with those in the Fig. 
2 of 35], we find that the value of fi m o constrained from the sample is the main reason of 
different evolutions of q(z) and Om(z). Because of the relative large Q m o obtained from 
Snla + BAO I data, the result is inconsistent with ACDM model at more than la level. 
This can be easily understood by expanding Om(z) for the CPL model at low redshift, 
Om(z) ~ 1 + u>o(l — n m o)- Even if we get wo close to —1, the large value of Q m o will make 
the CPL model inconsistent with ACDM model. However, if we impose a prior on f2 m o, for 
example, Q m o = 0.28 ± 0.04, then the result by fitting the Snla data or Snla+BAO I data 
to the CPL model, is consistent with the ACDM model. 



III. OBSERVATIONAL CONSTRAINTS ON PIECEWISE CONSTANT 
PARAMETRIZATION 

Although the CPL parametrization provides a useful tool to study the dynamical prop- 
erty of dark energy, the particular form of w(z) may impose a strong prior. Note that in 
a small enough redshift region, w(z) is approximately a constant, so we may divide the 
redshift into several bins, and parameterize w(z) as a constant in a particular redshift bin, 
this is the piecewise constant parametrization of w(z). If we have enough data, the red- 
shift range in each bin can be small enough, and the piecewise constant parametrization 
gives the true w(z). In other words, the piecewise constant parametrization is a model 
independent method. In practice, the number of redshift bins is finite, and the piecewise 
constant parametrization of w is an approximation of the true w(z), and it provides very 
useful information about the dynamical behavior of dark energy. In this section, we use 

6 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



FIG. 1: Reconstructed q{z) and Om(z) from Snla data, Snla+BAO II data and Snla+BAO 
II+CMB data using the CPL ansatz. The center solid lines are plotted with the best fit val- 
ues, and the shadows denote the la and 2a limits. The spatially flat ACDM model corresponds to 
a horizonal dashed line with Om(z) = 0.28 in the right panels. 

observational data to fit the piecewise constant parametrization of w(z). For the binning 



of the Constitution Snla data [39|], we apply the uniform, unbiased binning method [50 ]. 



We group the data into four bins so that the number of Snla in each bin times the width 
of each bin is around 30, i.e., NAz ~ 30, A" is the number of Snla in each bin, Az is the 
width of each bin. The choice of NAz ~ 30 for the Constitution Snla data results in four 
bins. The boundaries of the four bins are z\ = 0.22, zi = 0.53, z% = 0.85, z^ = 1.8 and 
z 5 extends beyond 1089. For the redshift in the range < z < z iy the equation of state 
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parameter is a constant, w(z) = W{. For convenience, we choose zq = 0. Due to the lack 
of observational data in the redshift range z — 1.8 — 1089, w(z) is largely unconstrained in 
this redshift range. For simplicity, we assume that w(z > 1.8) = —1. For a flat universe, 
the Friedmann equation becomes 



n 



e\z) = n m0 (i + z f + (i - n m0 )(i + z f 1+ ^ + Zi ^f^- Wi \ Zn . x < z <z n . (ii) 

i=l 

In this model, there are five free parameters fl m o, Wi, u>2, u>3 and u>4, let us denote them as 

= (Q m0 ,W 1 ,W2,W 3 ,W 4 ). 

In general, the equation of state parameters Wi in different bins are correlated and their 



errors depend upon each other. We follow Huterer and Cooray 221] to transform the covari- 
ance matrix of Wi to decorrelate the error estimate. Explicitly, the transformation is 

Wi^TTtjWj, (12) 



j 



where the transformation matrix T = V T A~ 1 ^ 2 V, the orthogonal matrix V diagonalizes 
the covariance matrix C of u>, and A is the diagonalized matrix of C. For a given i, Tij 
can be thought of as weights for each Wj in the transformation from Wi to Wj. We are 
free to rescale each Wi without changing the diagonality of the correlation matrix, so we 
then multiply both sides of the equation above by an amount such that the sum of the 
weights J? . is equal to one. This allows for easy interpretation of the weights as a kind 
of discretized window function. Now the transformation matrix element is Tij/^2, k Tik and 
the covariance matrix of the uncorrelated parameters is not the identity matrix. The i-th 
diagonal matrix element becomes {Y2jTij)~ 2 ■ In other words, the error of the uncorrelated 
parameters Wi is Oi = 

The likelihood for the parameters in the model and the nuisance parameters is computed 
using a Monte Carlo Markov Chain (MCMC). To observe the effect of the nuisance parameter 
Hq in the Snla data, we take two different approaches. In the first approach, we analytically 
marginalize Hq by using a flat prior 2^] • In the second approach, we take Hq = Hq as a free 
parameter in the MCMC code. The MCMC method randomly chooses values for the above 
parameters 0, evaluates y 2 and determines whether to accept or reject the set of parameters 
using the Metropolis-Hastings algorithm. The set of parameters that are accepted to the 
chain forms a new starting point for the next process, and the process is repeated for a 
sufficient number of steps until the required convergence is reached. Our MCMC code is 
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based on the publicly available package COSMOMC [52|. We give both marginalized and 
likelihood limits of the uncorrelated parameters Wj. The likelihood limit defines the region 
of parameter space enclosing a fraction / of the points with the highest likelihood as the 



3- 



The likelihood 



iV-dimensional confidence region, where / defines the confidence limit 
limit is very useful to assess the consistency with new data or theories. For comparison, we 
also fit the result by running the publicly available package WZBINNED [53] and the results 
are consistent with those by using the marginalized method. 

We first fit the parameters in the model by using the combined Snla+BAO II data, 
i.e., we calculate x 2 (#) = [A(0) - 0.469(0.96/0.98)~ - 35 ] 2 /0.017 2 + x 2 n (0,# o n ). By fitting 
the piecewise constant model to the data, we get x 2 = 459.2, the marginalized la estimate 
O m0 = 0.2851^011; an d the uncorrelated binned estimates of the equation of state parameters 
Wi are shown in Fig. |2j We see that the results marginalizing over Hq analytically are 
the same as those with H marginalized numerically. In other words, we can trust the 
result using the method of analytically marginalizing over Hq with a flat prior. From Fig. 
[21 it is clear that the likelihood limits are larger than the marginalized limits, especially 
for the parameters W3 and W4. If we take the likelihood limits, then the cosmological 
constant is consistent with the data even at la level. However, the cosmological constant 
is not consistent with the data at la level if we take the marginalized limits. Although 
the likelihood limits are useful information for the full MCMC sample, in the following, we 
quote the results obtained with marginalized limits. If we fit the flat ACDM model to the 
combined Snla+BAO II data, we get \ 2 = 466.4 and fi m0 = 0.286 ± 0.017. If we fit the 
flat CPL model to the data, we get \ 2 = 466.2, and the marginalized la constraints are 
Q m0 = 0.288±{j:oi6, w = -0.901^, and w a = -0.6j££. In Fig. El we show the marginalized 
probabilities of the parameters Wj. 

Now we add the BAO and WMAP5 data to the Snla data. The parameters 6 in the 
models are determined by minimizing x 2 (#> &bh 2 , h) = xL + Xbao + xlmbi nere h = ffi\/100- 



To use the BAO measurement (hereafter BAO III) from the SDSS data, we define 40|] 

2 fa n .a 7 \ f r s (z d )/D v (z = 0.2) -0.198 V f r s (z d )/D v (z = 0.35) - 0.1094 \ 2 
Xbao (0M b h,h)=^ J +^ J , 

(13) 



where the redshift z d is fitted with the formulae [54] 

1 + 0.659(fi m0 /i 2 ) a82T 



1291(fi m0 /i 2 ) ' 251 r , L2\fe 2 i riA\ 
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FIG. 2: The la and 2<r estimates of the four uncorrelated parameters Wj using the Snla+BAO II 
data. The parameter W4 is not well constrained. The error bars show la and 2a uncertainties with 
black solid lines and blue dotted lines, respectively. The top panels show the results by treating 
the nuisance parameter Hq as a free parameter, and the bottom panels show the results with Hq 
being analytically marginalized. The results in the left panels are the marginalized limits, while 
the results in the right panels are the likelihood limits. 



61 = 0.313{Q m0 h 2 )-° A19 [l + 0.607{Q m0 h 2 ) 0SU ] 
and the comoving sound horizon is 

dx 



b 2 = 0.238(ft m0 /i 



2\0.223 



r s (z) 



c s {x)E{xY 



(15) 



(16) 



where the sound speed c s (z) = l/y/3[l + R b /(1 + zj], and R b = 315000fi 6 /i 2 (2. 726/2. 7)~ 4 . 

To implement the WMAP5 data, we need to add three fitting parameters R, l a and z*, so 
Xcmfe(^) &bh 2 , h) = AxjCov _1 (xj, Xj)Axj, where Xi = (R, l a , z*) denote the three parameters 
for WMAP5 data, Axi = x% — x° bs and Cov(x{,Xj) is the covariance matrix for the three 
parameters [43|. The acoustic scale I a is 



Ia 



ird L (z* 



[1 + z*)r s (z*y 
10 



(17) 



_1 _o.9 -0.8 -0.7 

W 1 



-2 
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FIG. 3: The marginalized probabilities of the four equation of state parameters Wj. The dotted 
lines are the results using Snla+BAO II data only, the dashed lines are the results using the 
combined Snla, BAO III and WMAP5 data, and the solid lines are the results using Snla+BAO 
III+WMAP5+-ff (z) data. The nuisance parameter Hq in the Snla is analytically marginalized 
over. 



where the redshift z* is given by [55j 



2» = 1048[1 + 0.00124(n 6 /i 2 )- a738 ][l + gi (n m0 h 2 ) 92 } = 1090.04 ± 0.93, (18) 

_ 0.0783(^ 2 )-°- 238 0.560 
9l _ l + 39.5(rW) - 763 ' 92 ~ 1 + 21.1(JW) 1 - 81 ' ^ ' 

The shift parameter [43] 

f z * dz 

R(0, n b h 2 , h) = / — - = 1.710 ± 0.019, (20) 

Jo &\ z ) 
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When we add the BAO III and WMAP5 data, we add two more parameters VL^h 2 and the 
Hubble constant Hq = 100/i. Because the normalization of the luminosity distance-redshift 
relation is unknown, the nuisance parameter Hq in the Snla data is not the observed Hubble 
constant, and it is different from that in the BAO III and WMAP5 data. By fitting the 
data to the model, we find that the nuisance parameter Hq is around 65 km/s/Mpc in 
the Constitution data, while the Hubble constant Ho is around 72 km/s/Mpc. Therefore, 
we should treat the nuisance parameter Hq in the Snla data differently from the Hubble 
constant H in other observational data, and we should include both Hq and Hq {h) in the 
data fitting. We analytically marginalize over the nuisance parameter Hq in the Snla data 
as explained in 27] . If we treat the nuisance parameter Hq as the Hubble constant, we get 
X 2 = 465.9, and the uncorrelated estimates of Wj are shown in the top panels of Fig. HI If 
the nuisance parameter Hq in the Snla data is marginalized analytically, we get \ 2 = 462.3, 
n m0 = 0.283^q;q^> an d the uncorrelated estimates of Wj are shown in the bottom panels 
of Fig. HI By adding the BAO III and WMAP5 data, we add 4 more data points, and 
X 2 increases 3.1. With the help of the BAO III and WMAP5 data, the constraint on the 
equation of state parameters W; is improved, especially for W3 and W4. This point is also 
clear from the marginalized probabilities shown in Fig. [3j From Fig. HI we see that the 
results are different if the nuisance parameter Hq is treated differently. As explained above 



and in 



27] , we should treat the nuisance parameter Hq in the Snla data differently, so we 



analytically marginalize over Hq with a flat prior when fitting the Snla data. Again, we see 
that the marginalized limits give tighter constraints on the parameters than the likelihood 
limits do, so we quote the result in the bottom left panel as the fitting result. 

Finally, we add the data of the Hubble parameter H(z) at nine different redshifts from 
the differential ages of passively evolving galaxies obtained in 4f| and the three more recent 
data H(z = 0.24) = 79.69±2.32, H(z = 0.34) = 83.8±2.96, and H(z = 0.43) = 86.45±3.27 
by taking the BAO scale as a standard ruler in the radial direction 47j. To use these 12 
H(z) data, we define 

(8,h) = f2 [H ° bs{Zi) ; H{Zi)] \ (2i) 



12 

i=l ^hi 



where 07^ is the la uncertainty in the H(z) data. We also add the prior H = 74.2 ± 3.6 
km/s/Mpc determined from the observations with the Hubble Space Telescope by Riess et 



al. 



4§J. Now we have x 2 (#, V b h 2 , h) = X 2 n + xLo + xLb + xl- 
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FIG. 4: Same as Fig. except that we use Snla + BAO III + WMAP5 data here. 



If we treat the nuisance parameter Hq as the Hubble constant, we get x 2 — 500.5, and the 
uncorrelated estimates of Wj are shown in the top panels of Fig. If the nuisance parameter 
Hq in the Snla data is marginalized analytically, we get x 2 — 476.8, fl m o = 0.2711q;qqq, and 
the uncorrelated estimates of Wj are shown in the bottom panels of Fig. [5j By adding the 
H(z) data, we add 13 data points, and x 2 increases 14.5. As explained above, due to the 
difference between the nuisance parameter Hq in the Snla data and the Hubble constant 
H , we get much larger value of x 2 if we treat the nuisance parameter Hq in the Snla data 
as the Hubble constant, so we should use the results in the bottom left panel. If we fit the 
data to the fiat ACDM model, we get x 2 = 483.0 and Q mQ = 0.272 ± 0.011. If we fit the 



data to the fiat CPL model, we get x 2 = 482.5, Q m0 = 0.269±o;°o8, 



-0.97; 



-0.12 
-0.07' 



and 



w a = 0.031q75. These results suggest that the CPL model is consistent with ACDM model 
at la level, and ACDM model is consistent with the piecewise constant parametrization at 
2<T level. 

The marginalized probabilities of the parameters Wj are shown in Fig. [31 From Figs. [3] 
and El we see that the addition of the H(z) data improves the constraint on W4 by 70%. 
This suggests that the equation of state parameter at high redshift will be better constrained 
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with high quality data of H(z) in the future. 
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FIG. 5: Same as Fig. [21 except that we use the combined Snla, BAO III, WMAP5 and H{z) data 
here. 



IV. CONCLUSIONS 



Fitting the combined Snla and the BAO distance ratio (BAO I) data to the CPL model, 



-1.26 
-0.95 



we find that \ = 462.44, and the joint lcr constraints fl m o = 0.45_ 011 , wo = —0.13 
and w a = -12.2±\lj,. Fitting the combined Snla and the BAO A (BAO II) data to the 



CPL model, we find that \ 2 = 466.18, and the joint la constraints fl 



77(0 



0.291°;°!, w 



-0.901°;^ and w a 



-O.6I35. While the results obtained with the BAO I data are not 



consistent with the ACDM model at more than la level [35), the results obtained with 
the BAO II data are consistent with the ACDM model at la level (see Fig. [1]), and the 
constraints are tighter. So different BAO data give different results. The inconsistency lies 
mainly on the larger value of Q m Q obtained with the BAO I data. 

To constrain the property of dark energy using the observational data, we need to ap- 



ply model independent method. The piecewise constant parametrization of the equation 
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of state parameter w(z) of dark energy is somewhat model independent, we used the cur- 
rent observational data to study the property of dark energy with this model independent 
parametrization. Since the normalization of the luminosity distance-redshift relation is ar- 
bitrary, the nuisance parameter Hq in the Snla data is also arbitrary, and different from the 
observed Hubble constant Hq. We should treat it differently from the Hubble constant in 
other data, and we should include both Hq and Hq (h) in the data fitting. Otherwise, we 
may get wrong conclusions. If we treat the nuisance parameter Hq in the Snla data as the 
observed Hubble constant, then we may conclude that the flat ACDM model is incompat- 
ible with the combined Snla, BAO III and WMAP5 data or the combined Snla, BAO III, 
WMAP5 and H(z) data. However, by marginalizing over the nuisance parameter Hq ana- 
lytically in computing x 2 sn i the flat ACDM model is consistent with the current observation 
at 2a level. 

By fitting the combined Snla, BAO III, WMAP5 and H(z) data to the flat ACDM model, 
we get x 2 = 483.0 and fi m0 = 0.272 ±0.011. If we fit the data to the flat CPL model, we get 
X 2 = 482.5, and the marginalized la constraints f2 m0 = 0. 2691^08, = ~~ 0-97lao7' anc ^ 
w a = 0. 03^0^5. If we fit the model with piecewise constant parametrization to the data, we 
get x 2 = 476.8, and the marginalized la constraints Q m0 = 0.27llo;°o6, W i = _0 - 93 ± °- 06 > 
m = -1-131°;}°, m = -0.1218;|, and W 4 = -0.85!?;}jj. 

By adding the H(z) data to the Snla, BAO III and WMAP5 data, the constraint on 
the parameters is greatly improved, especially for the equation of state parameter at high 
redshift where the number of Snla data is small. The result suggests that the equation of 
state parameter at high redshift will be better constrained with high quality data of H(z) 
in the future. 
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